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Abstract 

The temperature phase transition in the A^-component scalar field the¬ 
ory with spontaneous symmetry breaking is investigated using the method 
combining the second Legendre transform and with the consideration of 
gap equations in the extrema of the free energy. After resummation of all 
super daisy graphs an effective expansion parameter, (l/ 2 A^)^e^ appears 
near Tc for large N. The perturbation theory in this parameter account¬ 
ing consistently for the graphs beyond the super daisies is developed. A 
certain class of such diagrams dominant in 1/N is calculated perturba- 
tively. Corrections to the characteristics of the phase transition due to 
these contributions are obtained and turn out to be next-to-leading order 
as compared to the values derived on the super daisy level and do not alter 
the type of the phase transition which is weakly first-order. In the limit N 
goes to infinity the phase transition becomes second order. A comparison 
with other approaches is done. 


1. Introduction 

Investigations of the temperature phase transition in the iV-component scalar 
held theory (O(A^)-model) have a long history and were carried out by either 
perturbative or non perturbative methods. This model enters as an important 
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part unified field theories and serves to supply masses to fermion and gauge fields 
via the mechanism of the spontaneous symmetry breaking. A general believe 
about the type of the symmetry restoration phase transition is that it is of sec¬ 
ond order for any values of N (see, for instance, the text books |ll-0). This 
conclusion results mainly from non perturbative analytic and numeric methods 
[il”lll' opposite, a hrst order phase transition was observed in most pertur¬ 
bative approaches [P]-|]l4|. An analysis of the sources of this discrepancy was 
done in different places, in particular, in our previous papers |T3-|T3 devoted 


to the investigation of the phase transition in the (9(A^)-model in perturbation 
theory (PT). Therein a new method has been developed which combines the 
second Legendre transform with consideration of the gap equations in the ex¬ 
trema of the free energy functional. This allows for considerable simplihcation of 
calculations and for analytic results. The main of them is the discovery in the 
so-called super daisy approximation (SDA) of an effective expansion parameter 
e = near the phase transition temperature T^. All quantities (the particle 
masses, the temperatures T+,T_) are expandable in this parameter. The phase 
transition was found to be weakly hrst order converting into a second order one 
in the limit A^ —^ cx). The existence of this small parameter improves the status 
of the resummed perturbative approach making it as reliable as any perturba¬ 
tive calculation in quantum held theory. For comparison we note that formerly 
there had been two problems with the perturbative calculations near the phase 
transition. First, unwanted imaginary parts had been observed. Using functional 
methods we showed in |^, that they disappear after re-summing the super 
daisy graphs. The second problem was that near the masses become small 
(although not zero) compensating the smallness of the coupling constant hence 
making the ehective expansion parameter of order one. 

In the present paper we construct the PT in the effective expansion parameter 
^ 0(A^)-model at large N. It uses as input parameters 

the values obtained in the SDA. As an application we calculate corrections to 
the characteristics of the phase transition near T^. which follow from taking into 
account all BSDA graphs in order T Since the masses of particles calculated in 
the SDA are temperature dependent, we consider in detail the renormalization at 
hnite temperature of the graphs investigated. It will be shown that the counter 
terms renormalizing these graphs at zero temperature are sufficient to carry out 
renormalizations at hnite temperature. 

The paper is organized as follows. In sect.2 we adduce the necessary in¬ 
formation on the second Legendre transform and formulate the BSDA PT at 
temperatures T ~ T^. In the next section we calculate the contribution to the 
free energy of the graphs called the “bubble chains” having a next-next-to-leading 
order in A. In sect.4 the renormalization is discussed. The corrections to the 
masses and other parameters near T^. are calculated in sect.5. The last section is 
devoted to discussion. 
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2. Perturbation theory beyond super daisy ap¬ 
proximation 


In this section we develop the PT in the effective expansion parameter e = -^73 
for the graphs BSDA in the frameworks of the second Legendre transform. Con¬ 
sideration of this problem is quite general and independent of the specific form 
of the Lagrangian. So, it will be carried out in condensed notations of Refs. 


I1-0 


The second Legendre transform is introduced by representing the connected 
Green functions in the form 


W = R[0] + ^Trlog/3 - + W^, 


( 1 ) 


where S'[0] is the tree level action, (3 is the full propagator of the scalar N- 
component field, A is the free field propagator, W 2 is the contribution of two 
particle irreducible (2PI) graphs taken out of the connected Green functions and 
having the functions j3{p) on the lines. The symbol ”Tr” means summation over 
discrete Matsubara frequencies and integration over a three momentum (see for 
0 - 0 ). 


more details Refs. 


The propagator is related to the mass operator by the Schwinger-Dyson equa¬ 


tion 


r'(p) = A-i(p)-E(p), 
5W2 


E(p) = 2 


6f3{p) • 


( 2 ) 

(3) 


The general expressions (|^) and (^) will be the starting point for the construction 
of the BSDA PT. Galculations in SDA have been carried out already in |]I^ - 


1 ^ and delivered the masses of the fields and the characteristics of the phase 


transition: Tc - transition temperature, and T+,T_ - upper and lower spinodal 
temperatures. These parameters will be used in the new PT as the input parame¬ 
ters (zeroth approximation). Gontributions of BSDA diagrams will be calculated 
perturbatively. 

First let us write the propagator in the form 


/?(p) = /?o(p) +/?'(p), 


(4) 


where /do is derived in the SDA and f3' is a correction which has to be calculated 
in the BSDA PT in the small parameter e = for large N. The 2PI part can 
be presented as 


W 2 = WsD + w' = WsdIPo + 0] + W!,[l3o + 0] (5) 
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and assuming (3' to be small of order e value we write 


WsD = WsDm + ^-^^^(^' + 0{e^), ( 6 ) 

dPo 

Wi = + + (7) 

In the above formulas the squared brackets denote a functional dependence on 
the propagator. The curly brackets as usual mark a parameter dependence. In 
such a way other functions can be expanded. For we have 


r‘ = 


( 8 ) 

( 9 ) 


where S'[/3o] = jg ^ correction following from the 2PI graphs, So[/3o](p) 


is the super daisy mass operator. In a high temperature limit within the ansatz 
adopted in Refs.|[T^-[p^ + M^) it looks as follows 


P \p)=p‘^ + Ml - S'[/3o](p), 


( 10 ) 


where Ml is the held mass calculated in the SDA as the solution of the gap 
equations. 

In a similar way, the free energy functional can be presented as 

W = + W' (11) 


with 


and 


+ ^TrlogPo - ^Tr/^oA'^ + WsoW 


QiK-l , ^^SD[Po\ni , Tj/r/3 1 , q! a-l , 2 




7('+lV2|/3„] + -rr/3'/3„-' + 0(€"). 


Taking into account that P' =2 one hnds 

lrrW„-> = 

and hence 


W' = W'^[Po]. 

Thus, within the representation (^) we obtain for the IF functional 


( 12 ) 

( 13 ) 

( 14 ) 

( 15 ) 

( 16 ) 
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where is the expression (|l]) containing as the Ph 2 [/ 9 o] the SDA part only 

and the particle masses have to be calculated in the SDA. The term 1^2 [/^o] 
corresponds to the 2PI graphs taken with the /do propagators on lines. 

From the above consideration it follows that perturbative calculations in the 
parameter e derived in the SDA within the second Legendre transform can be 
implemented in a simple procedure including as the input masses of propagators 
/d the ones obtained in the SDA. Different types of the BSDA diagrams exist. 
They can be classihed as the sets of diagrams having the same orders in A from 
the number of components. So, it is reasonable to account for the contributions 
of the each class by summing up all diagrams with the corresponding specihc 
power of 4. A particular example of such type calculations will be done below. 


3. Expansion near Tc 


In this section we shall calculate a hrst correction in the effective expansion 
parameter e = -^^^73 for large N ai T ~ Tc. 

Before to elaborate that, let us take into consideration the main results on 
the SDA which have to be used as a PT input. In Ref. |T3 it was shown that the 


masses of the Higgs and the Goldstone helds near the phase transition 
temperature are ( for large N) 


(17) 


A/flO) _ 

AT+ 

^ 1 

1 + ^ 

IVlri 

dvr 

42Ar)V3 

2N ■■'A 


AT+ 

( 1 



27r 

l(2Ar)2/3 

2N ^ 


where A is a coupling constant, upper script zero means that in what follows these 
masses will be chosen as zero approximation, m - initial mass in the Lagrangian. 
The upper spinodal temperature T+ is close to the transition temperature Tc ^ ^ 
. We adduced here the masses for the T+ case because that delivers simple analytic 
expressions. Results for other temperatures in between, T_ < T < T+, are too 
large to be presented here. Note also that the temperatures T+,T_ in the SDA 
are related as (see Refs. 


13 , \14 


— = 1 + - - - 

T_ 167r2 (2Ar)2/3 


+ ..., 


(18) 


and T_ = y x{n+ 2 ) "^' ff^^^e it is clear that the transition is of a weakly hrst-order 
transforming into a second order one in the limit N ^ 00 . 

With these parameters taken into account an arbitrary graph beyond the SDA 
having n vertexes can be presented as 




3C’-2Lt/ _ \1/ 
n 




( 19 ) 
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Here the notation is introduced: C = L- n + l- number of loops, L - number 
of lines. Since we are interested in diagrams with closed loops only, the relation 
holds: 2n = L. The vertex factor Vn comes from combinatorics. The right-hand- 
side of Eq.(|T^) follows when one shifts the three dimensional momentum of each 
loop, p —> Mp, and substitutes instead of M the mass Eq. fpTf) to have the 
’worst case’ to consider. In this estimate the static modes (Z = 0 Matsubara 
frequency) of propagators were accounted for. One may wonder is it sufficient at 
Tc? The positive answer immediately follows if one takes into consideration that 
side by side with the three-momentum rescaling the temperature component of 
the propagator is also shifted as T —^ ~ Hence it is clear that at N 

goes to infinity a high temperature expansion is applicable and the static mode 
limit is reasonable. 

As it follows from Eq. (p!9D , A is not a good expansion parameter near Tc 
whereas is the one because it enters in the power of the number of vertexes 
of the graph. Of course, we have to consider the graphs beyond the SDA. That 
is, all diagrams with closed loops of one line (tadpoles) have to be excluded 
because they were summed up completely already at the SDA level. Note also 
the important factor ^ coming from the rescaling of three-momentum with the 


mass M, 


( 0 ) 


0 


Before to proceed we note the most important advantages of resummations 
in the SDA [|^|,|H]: 1) There is no an imaginary part in the extrema of the 


free energy. 2) The simple ansatz for the two-particle-irreducible Green function 
= p^ + is exact in this case. Here, p is a four-momentum, M is a mass 
parameter which is determined from the solution of gap equations. 3) As it is 
known for many years, is well determined by this approximation and it is 
not altered when further resummations are achieved. 4) The existence of the 
expansion parameter e = near Tc. 

Since the estimate (|l^) assumes the rescaling of momenta p Mp the same 
procedure should be fulfilled in perturbation calculations in the parameter e. 
They are carried out in the following way. First, since A is not the expansion 
parameter it can be skipped. Second, only the BSDA diagrams have to be taken 
into consideration. Moreover, since is a large number it is convenient to sum up 
series having different powers in A. Third, rescaling p —> Mp can be done before 
actual calculations. In this case one starts with the expressions like in Eq.([I^) 
(for diagrams of 0-sort to consider). When p-particles are included one has to 
account for the corresponding mass value and rescale the momentum accordingly. 
Fourth, the temperature is fixed to be Tc. In other words, one has to start with 
the expressions like in Eq. (0). 

To demonstrate the procedure, let us calculate a series of graphs giving leading 
in A contribution F' to the free energy F. Then, the total to be T = -|- T\ 

where is the result of the SDA. 

The ’bubble chains’ of 0-field are the most divergent in N and other diagrams 


6 










have at least one power of N less. So, below we discuss the 0- bubble chains, 
only. The contribution of these sequences with the mass and various n is 
given by the series 

\ 1 CO -I T A 

o* = of + ^ y; _^r.,( 6 EW(p)iv^/»)" (20) 


Here, is the contribution of the ’’basketball” diagram, is the diagram 

of the type 

+ p) (21) 

and the power of the parameter A is written explicitly. 

Now, let us sum up this series for large hxed N. The vertex combinatorial 
factor for the diagram with n circles is 0 


K = + ilr-\N - 2)]. (22) 

The leading in N term in the 14i is ~ {N + 1/3)"’ and we have for the series in 

Eq.dl) 

\ 1 CO -| 

JH H + l/W)"]. (23) 

n=3 ^ 

To sum up over n we add and subtract two terms: —^+ 1/A^), and 
il-Ej>(p)A'^''“(A' + l/A')]^. We get 


-^TT,llog(l + T.ll\p)N^/^N+l/N)) (24) 

- + l/A') + + 1/A'))"]- 

To hnd all together we insert for the hrst term 


D 


m i u (iv"-i) r, <-4bp))" 

4 iV2 ^ iV2/3 


Then the limit N goes to inhnity has to be calculated. We obtain hnally: 

A' 


= -^rr>j(l + E«(p)lV"l’) 

^Tr,Ej>(p) - j^Tr,(E5.‘’(p))= 


(25) 


(26) 


Thus, after summing over n the limit N goes to inhnity exists and the series is well 
convergent. Eq.(PBD gives the leading contribution to the free energy calculated 
in the BSDA in the limit N goes to inhnity. This is the hnal result, if we are 
interested in the leading in jr correction. 
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To obtain the squared mass, M^, one has to sum up the value calculated in 

the SDA in accordance with a general expression 

for the mass []^. By substituting the propagators one is able to hnd leading in 
corrections to free energy and masses at temperatures close to Tc. In this 
way other quantities can be computed. 

The most important observation following from this example is that the limit 
N goes to inhnity does not commute with summing up of inhnite series in n at 
n —>■ oo. It is seen that the hrst two terms in the Eq.([^) can be neglected as 
compared to the last one which is leading in 1/N. As it is occured, this contri¬ 
bution is of order that is smaller then the value of ~ A^“A3 

determined in the SDA. Note that the correction is positive and the complete 
held mass is slightly increased. We shall calculate the value of the mass in sect.5 
simultaneously with other parameters of the phase transition. Contributions of 
other next-next-to-leading classes of BSDA diagrams can also be obtained per- 
turbatively. 


4. Renormalization of vortexes and bubble chains 


Calculations carried out in the previous section deal with the unrenormalized 
functions. But the question may arise: whether the graphs, which are series of 
the function with the temperature dependent mass are renor¬ 

malized by the temperature independent counter terms as one expects at hnite 
temperature? We shall consider in detail this question for the non symmetric ver¬ 
texes and the graphs of ’’bubble chain circles ” with n M^) 

insertions. These are of interest for computations in the previous section. We will 
show that this is the case for both of these objects. They contain, correspond¬ 
ingly, n — 1 and n insertions of (p, M^), p = pa + Pb a momentum incoming 
in the one-loop vertex M<^). The functions are calculated from 

^absdip) by means of contracting the indexes and integrating over internal mo¬ 
mentum to form one extra M^) term. Note that the Goldstone held bubble 
chains give a leading in A contribution among the BSDA diagrams. 

To carry out actual calculations we adopt the 0(V)-model with the notations 
introduced in Refs. 0, 0- In the restored phase, we have N scalar helds with 
the same masses Mr{T). In the broken phase, there is one Higgs held p with the 
mass Mrj and N — 1 Goldstone helds 0 having the mass (0 derived in the 
SDA. 

Let us hrst consider the one-loop vertex of the 0-held in the s-channel, s = 
= {Pa + PbY, CL, b mark incoming momenta. 


^abcidi 2 2 T r^abc\d^ i 


U 


(27) 
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where the notation is introduced: p = - expansion parameter in the 0{N)- 

model, Cl = — |) is a combinatorial factor appearing at the symmetric 

tensor ^^ab<^cidiy ^abcidi ^(^ab^Cidi -|- <^aci^bdi “ 1 “ ^adi^bci') is the tree 

vertex in the 0(iV)-model. Subscript 1 in the indexes ci,di counts the number 
of loops in the diagram. The diagram with n loops ( insertions) has the 

form 


■^(n+l) 

ClbCndn 


p{n+l) 

2 


iY^iCnSabc 


nd 


71 


_L 1 

^ V 2 / ^ abcndji / ’ 


(28) 


where now Cn = and other notations are obvious. 

At zero temperature, the function has a logarithmic divergence which 

in a dimensional regularization exhibits itself as a pole ~ -f, e = d —4. We denote 
the divergent part of S^(p) as Di. To eliminate this part we introduce into the 
Lagrangian the counter term C 2 of order p^\ 


C 2 


2-^1 (2) 
P 2 


(29) 


where 


(Cis + ^V^)abcidi- Thus, the renormalized one-loop vertex is 


y{2) ^ ^ 2 ^ren. (2) 

^ abcidi ryr r\ ^abcidi 


(30) 


with - Di. 

In order three diagrams contribute. The hrst comes from the tree vertexes 
and contains two insertions. Two other graphs are generated by 
and the counter term vertex 6 * 2 . Each of them has one insertion. The sum of 
these three diagrams is 


= 

abc2d2 




(31) 


The terms with the products Di cancel in the total. To have a hnite expres¬ 
sion one has to introduce a new counter term vertex 


C 3 


rP^DtP 


2 ., (3) 


abc2d2 • 


(32) 


It cancels the last independent of divergence. Thus the renormalized vertex 
y(3) jg given by the hrst term in the expression (pl|). 

This procedure can be easily continued with the result that in the order pC+^) 
one has to introduce the counter term of the form 


C-n+l P 


n+1 



) ^abcndn ■ 


(33) 


The hnite vertex calculated with all types of diagrams of the order p^~^^ looks as 
follows 


1 

^n+l _ _ n +1 ^ C^re^\n (n+l) 

^ abcndn 2 " '^2 ^ahcndn' 


(34) 
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Above we considered the s-channel diagram contributions to the vertex . This 


is sufficient to study the leading in A correction (^ 6 [) of interest. To have a 


symmetric renormalized vertex one has to add the contributions of the t- and 
M-channels and multiply the total by 

Now it is easy to show that the counter terms Cn rendering the hniteness of 
the vertexes at zero temperature are sufficient to renormalize them at hnite tem¬ 
perature. Really, M(T),T) can be divided in two parts, M(T),T = 
0) = M(T)) corresponding to held theory and M(T),T) - the statis¬ 

tical part. The divergent part Di and the counter terms Cn are independent of 
mass parameters. So, to obtain the renormalized vertex at hnite temperature 
^abcndniP^ T) one has to sum up the same series of the usual and the counter term 
diagrams as at zero temperature and substitute M(T), T) = M{T)) -|- 

S;^(p, M{T),T). Again, the terms of the form [S^(p, are canceled 

and we hnd the hnite expression 

KE.y(p.r,M(r)) = ^ + (35) 


For this procedure to hold it is important that Di is logarithmically divergent 
and independent of mass. So, the held theoretical part of T, M(T)) 

as well as the statistical part is renormalized by the same counter terms as the 
vertex at zero temperature. In the BSDA PT we have to use the mass 

Let us turn to the functions They can be obtained from the vertexes 

1/T, M(T)) in the following way. One has to contract the initial and the 
hnal indexes a, b and c„, dn and form two propagators Po which after integration 
over the internal line give the term The combinatorial factor of this diagram 
is ^- 7 . 

n+1 

We proceed with the function at zero temperature. In the order two 
diagrams contribute. One includes three ordinary vertexes 


s?hp)= 2^v'z,{^y, 


(36) 


and two other diagrams containing one vertex and the counterterm vertex 

^^2 \P) qP ^abcd^cdaby o / 9 * 


(3) 


(37) 

_ n3 — -(3) _ 

— dj, — — 


The contraction of indexes in Eq.(|37|) gives the factor v 

+ (|)^(^ ~ 2 )]. Again, substituting -|- Di one can see 

that the terms with products are canceled in the sum of expres¬ 

sions (^61) and (^71) . To obtain a hnite we introduce the counter term of order 


C(3) = lp^£l 
3^2 


n3 


( 38 ) 
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( 39 ) 


After that the renormalized circle is 

slfjp) = 

Note that since does not depend on any parameter, it can be omitted as well 
as the divergent term in the expression In other words, there is no need to 
introduce new counter terms into the Lagrangian in order to have a hnite 
and the counter term renormalizing the vertex are sufficient. 

This procedure can be continued step by step for diagrams with any number 
of insertions. The renormalized graph looks as follows, 

s'SKp) = (“) 

where the factor coming from the contraction of the 


d: = 


N -\- 1 r/iV + l\n-l 




(41) 


Now, it is a simple task to show that the counter terms renormalizing vertexes 
and circles 5^”^ at zero temperature are sufficient to obtain hnite S^'^\p,T) 
when the temperature is switched on. This is based on the property that the 
temperature dependent graph S^(p, M, T) can be presented as the sum of the zero 
temperature part and the statistical part independently of the mass term entering. 
Then it is easy to check that the divergent terms of the form [S^(p, M, T = 
0) + M, T)]*(Zli)"“^ are canceled when all the diagrams of order forming 
the circle Sl^^ {p,T) are summed up. 

Thus we have shown, the renormalization counter terms of leading in 1/iV 
BSDA graphs calculated at zero temperature being independent of the mass 
parameter entering renormalize the T) functions at hnite temperature. 

This gives the possibility to construct a PT based on the solutions of the gap 
equations. In fact, just the series of functions are of interest at the transition 
temperature which has to be considered as a hxed given number. Other 
parameters of the phase transitions can be found by means of some iteration 
procedure of the investigated already gap equations. 


5. Corrections to the parameters of the phase 
transition 

Having obtained the leading correction to the free energy (PB|) one is able to hnd 
perturbatively the characteristics of the phase transition near Tc. We shall do 
that for the limit N ^ oo. 
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First let us calculate the corrections to the Higgs boson mass, AM^j, and the 
Goldstone boson mass, AM^, due to term The starting point of this 

calculations is the system of gap equations derived in Ref. (Eqs. (28), (30) 
and (45)) and Ref. [jT^ (Eq. (28)). Here we write that in the form when only 


the term containing is included. 


m2 

7] 

2 


= m 


A 


(0) 


A 


iV-1, 

N ' 




40) 


6D, 


(42) 


-V ^ N- 16 p^{ 0 )' 

Remind that these equations give the spectrum of mass in the extrema of the free 
energy in the phase with broken symmetry. The complete system (Eqs. (45) in 
Ref. [T^) contains other terms which have higher orders in ^ and were omitted. 

Here m is the mass parameter in the Lagrangian. The functions , E® are 
the tadpole graphs with the full propagators Pn,P(t> on the lines. For the ansatz 


adopted in Refs.ll^, 
asymptotic expansion 


, they have at high temperature the 


= Tr/3,^* = 


M„,4,T 


+ 


12 47r 

where dots mark next-next-to-leading terms. Within Eqs.(^) - 


(43) 

(without the 


term) the masses o have been derived in the limit of large N. 

Now we compute the last term in the Eq. (^) for large N. In this case the last 
term of in the Eq.(^6D is dominant and calculating the functional derivative 
we hnd 

r(0)^ _ ^ ^ 1 


= 




(44) 


N-1 6(3^{0) N-1 

The sunset diagram entering the right-hand-side can be easily computed to give 


13 


3M 


(0) 




Thus for /(M^°^) we obtain in the large N limit 

A T2 


N^/3 32tA ^'",r(2Af)2/4’ 


(45) 


(46) 


where the mass M® (|^ was inserted and T = has to be substituted. Here 
we again turn to the T+ case to display analytic results. 

Since /(M®) is small, it can be treated perturbatively when the masses 
and are calculated. Let as write them as 

(47) 


Mrj = M^+x, 

= Mf+1/ 
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assuming x, y to be small. Substituting these in the equations 
preserving linear in x, y terms we obtain the system 


2Mf 1/ = 


AT 


2.N 

3AT \{N - l)T 
X + T, y 


= 2,^ 


and 


(48) 


2nN 


of linear equations. Its solutions for large N are 


1 T 


X = 


y = 


+ 


1 


(l-2ln 


3 327r Ar2/3 

1 2^/3 T+ 1 / 

2 ^ m iVV^ 


3^^ N 
7r(2iV)2/3J’ 

3^^ 

7r(2Ar)2/3 


(49) 


2 In 


As one can see, these corrections are positive numbers smaller than the masses 
(0 calculated in the SDA, as it should be in a consistent PT. Notice that the 
correction x is larger as compared to the next to leading term in the SDA (l/iV). 
So, the BSDA graphs deliver the next-to-leading correction. The value y is of the 
same order as the next-to-leading term in the SDA. 

In a similar way the correction to the mass in the restored phase, AMr, can be 
calculated. In this symmetric case the all components have equal masses which 


in the SDA within the representation (^3[) are the solutions of the gap equation 
(Eq. (36) in Ref.|]I^ and Eq.(27) in Ref. H ) 


= 


2 A(iV + 2)T 


12iV 


A(iV + 2)T 


SttN 


(50) 


It has a simple analytic solution 


,„,_ A(iv + 2)r M(iv + 2)r .a_ 

StiN V V Htt/V J 


SnN 


X{N + 2)T^ 
12iV ’ 


(51) 


where again the upper script zero means the SDA result. The contribution of the 
bubble graphs (^) corresponds either to the broken or to the restored phases. 
The only difference is in the number of contributing held components. In the 
broken phase this is iV — 1 and in the restored it is N. But it does not matter for 


large N. So, the function which has to be substituted into Eq. (^) is the one in 


Eq. (1^) with the replacements N — 1 ^ N and —>• . 

To calculate the mass correction z which is assumed to be small we put Mr = 
-|- 2 ; into Eq.(p® and hnd in the limit N ^ 00 : 


z = 


/(M<°>) 

2M)“' + y 

47r 


(52) 
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We see that the correction has a small positive value. By requiring that Mr to be 
positive the temperature T_ can be simply calculated. It is clear that T_ < 
and BSDA graphs decrease slightly the lower spinodal temperature. 

Now we calculate the correction to T+. To do that let us turn again to the 
system (^) with the expressions (|^), (^61) been substituted. From the second 
equation of the system we hud 


Mr] — Md) + 


2'kN 

'W 


(MJ - /(Mf)) 


and insert this into the hrst equation to have 

M^) = 


C-^)\K - 2Mlf) + Mj + ^^(Ml 


2m^ 


XT 

A(iV + 2)T2 3AT/ ^ 27rA^...2 


(53) 


(54) 


6N 


2'7iN 


XT 


/)) + 


A(iv- i)r 

2tiN 




where the linear in / terms are retained. This fourth order algebraic equation 
can be rewritten in the dimensionless variables /i = and g = (^^)^/ as 

follows: 


F(/i) = 0, (55) 

F(/i) = + 2fi^ — 2fi'^ — h — {N + 2)fi — 2g{g.'^ + fi). 


Here the /i-independent function is 


h = — 3g = {2m? 


X{N + 2)T‘^ ^^2'kN 


6N 


XT 


^9 


(56) 


and as before the SDA part is marked by the upper script zero. Remind that this 
equation determines the Goldstone field masses in the extrema of the free energy. 
It has two real solutions corresponding to a minimum and a maximum. The 
equation (|55[) has roots for any N and T which are too large to be displayed here. 
To have simple analytic results and investigate the limit of large N we proceed 
as in Ref.|T^ and consider the condition for the upper spinodal temperature. In 
the case of T = T+ the two solutions merge and we have a second equation, 


F (/i) — /^^ + ~ 2^ ~ 2^ ~ 2^ ~ 

It can be easily solved for large N: 

Here the hrst two terms of asymptotic expansion are included that is sufficient 
for what follows. Notice that since ^f-dependent terms are of next-next-to-leading 
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order they do not affect the main contributions. The function h can be calculated 
from Eq.f^Sl) and its first two terms are 




/ \ 4/3 

-3(t) 


7 /iV \ 2/3 
2M / 


+ ... 


The temperature T+ computed from Eq. 


can be rewritten in the form 


- + 


^27rA^/27r^(A^ + 2)A^ ,\ 

= ( 3A 


- 1/2 


(59) 


(60) 


where the values of Eq.(|5^ ) and g Eq. (^) have to be substituted. Again, 
since g is a small next-next-to-leading term the solution can be obtained pertur- 
batively. We find 


T^ = T 


(0) 


0 + 


9A (1 


9 Jn 3\/3A ' 

^ '''''7r(2Af)2/3, 


327r2 


Ar5/3 


(61) 


where the value = T^\l -|- must be inserted. From this result 

it follows that the upper spinodal temperature is slightly increased due to the 
BSDA contributions. But this is next-next-to-leading correction to the of 
the SDA. 

Thus, we have calculated the main BSDA corrections to the particle masses 
and the upper and lower spinodal temperatures in the PT. We found that 

r_ is decreased and T+ is increased as compared to the SDA results due to the 
leading in T BSDA graphs - bubble chains. So, the strength of the first-order 
phase transition is slightly increased when this contribution is accounted for. 
However, this is a next-to-leading effect. In such a way we prove the results on 
the type of the phase transition obtained already in the SDA [O, [IT3 . 


6. Discussion 


As the carried out calculations showed, the phase transition in the 0(A^)-model at 
large finite N is weakly of first-order. It becomes a second order one in the limit 
N ^ oo. This conclusion has been obtained in the SDA and was proved in the 
perturbative calculations in the consistent BSDA PT in the effective expansion 
parameter e = (;^)^^^- This parameter appears near the phase transition temper¬ 
ature Tc at the SDA level. Let us summarize the main steps of the computation 
procedure applied and the approximations used to derive that. 

In Refs.||12|| - as a new method the combination of the second Legendre 
transform with considering of gap equations in the extrema of a free energy func¬ 
tional was proposed. This has simplified calculations considerably and resulted in 
transparent formulas for many interesting quantities. Within this approach the 
phase transitions in the 0(1) - and 0(A^)-models were investigated in the super 
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daisy and beyond approximations. To have analytic expressions a high tempera- 
tnre expansion was systematically applied and the ansatz for the fnll propagators 
has been nsed. This ansatz is exact for the SDA which snms 
np completely the tadpole graphs. Jnst within these assnmptions a hrst order 
phase transition was observed and the effective parameter e has been fonnd. In 
terms of it all interesting characteristics can be expressed in the limit of large 
N and pertnrbation theory at T ~ constrncted. This solves an old problem 
on the choice of a zero approximation for pertnrbative compntations near T^. (if 
snch exist). We have shown explicitly here that these are the SDA parameters 
taken as the inpnt values for BSDA resummations. Clearly, this does not change 
qualitatively the results obtained in the SDA. 

Two important questions should be answered in connection with the results 
presented. The hrst is on the relation with other investigations where a second 
order phase transition was observed (see, for example, the well known literature 
i-i). and in fact this is a general opinion. The second is the non zero mass 
(P!7D of the Goldstone excitations in the broken phase at hnite temperature that 
was determined in SDA by solving of gap equations ||I3[ 


What concerns the hrst question, we are able to analyze the papers [Q, p [|T^ , 


1(;],|17|, where analytic calculations have been reported. Partially that was done 


in Refs. [|T^ , |l^ . We repeat that here for completeness. First of all we notice that 
there are no discrepancies for the limit N ^ oo where all the methods determine 
a second order phase transition. The discrepancy is for large hnite N. 

In Ref. ||15|| the renormalization group method at hnite temperature has been 
used and a second order phase transition was observed. Results obtained in 
this approach are difficult to compare with that found in case of a standard 
renormalization at zero temperature. This is because a renormalization at hnite 
temperature replaces some resummations of diagrams which remain unspecihed 
basically. In Refs.J^, a non-perturbative method of calculation of the ef¬ 
fective potential at hnite temperature - an auxiliary held method - has been 
developed and a second order phase transition was observed in both, the one- 
and the A^-component models. This approach seems to us not self-consistent be¬ 
cause it delivers an imaginary part to the ehective potential in its minima. This 
important point is crucial for any calculation procedure as a whole. Really, the 
minima of an ehective potential describe physical vacua of a system. An imagi¬ 
nary part is signaling either the false vacuum or the inconsistency of a calculation 
procedure used. This is well known beginning from the pioneer work by Dolan 
and Jackiw [|T^ who noted the necessity of resummations in order to have a real 


ehective potential. In our method of calculation this requirement is automatically 
satished when the SDA diagrams are resummed in the extrema of free energy. 
This consistent approximation is widely used and discussed in diherent aspects 
in literature nowadays P^-[pi[ . 

Let us discuss in more detail the results of Refs.Q], ^ where an interesting 
method - the method of average potential - was developed and a second order 
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phase transition has been observed for any value of N. To be as transparent 
as possible we consider the equation for the effective potential derived in Ref.p 
(equation (3.13)) 

U’ (p, T) = X{a^ + p-b ^JU'{p,T )), (62) 

where the notation is introduced: U'{p,T) is the derivative with respect to p of 
U{p,T), p = is the condensate value of the scalar held, A is the coupling 

constant. The parameters and b are: 


b 


(T^ - 
NT 



Sir 


24:771% 

A^’ 


(63) 


Here rriR and A/j are renormalized values, and Xr = A(1 + (see 

for details the cited paper). 

Considering the temperatures T ~ T^r and p « the authors have ne¬ 
glected U'{p,T) as compared to the {p, T) in Eq. (|6^ and obtained after 
integration over p (formula (3.16) of Ref. P) 


^(P, T) 


7 r 2 .T2 - T2 2 1 87^2 T2 - Tl 2 1 1 3 

9 ( T ) ^ /V 3 P ^ 3 N ' 


(64) 


Near T ~ T^r this effective potential includes the p^{4>^) term and predicts a 
second order phase transition. That is the main conclusion. 

However, this analysis is insufficient to distinguish between a second and a 
weakly-hrst-order phase transitions. Actually, the temperature Ter determined 
from the initial condition IT {0, Ter) = 0 in the former case corresponds to the 
temperature T_ in the latter one. To determine the type of the phase transition 
one has to verify whether or not a maximum at hnite distances from the origin in 
the p-plane exists. To check that the equation (|62|) has to be integrated exactly, 
without truncations. This is not a difficult problem if one hrst solves Eq. ([6^) 
with respect to U {p,T), 


U {p,T) — X(a'^ + + p — b ^A(a2 + p) -\ — —'j , 


and then integrates to obtain 


TT7 _ f^7rr2 , i 2 

U(p,T) - (A(T -T„)-+^^)p + -p 


2XTN 

3 Stt 


{{r 


N 


2^- AT2iV2 , ^3/2 

"’'^24 2567r2 


A 3/2 

+ P) ) 


(65) 


( 66 ) 


where the values of a and b Eq.(|6^) are accounted for. We see that in contrast 
to Eq.(|6^ this potential includes a cubic term which is responsible for the ap¬ 
pearance of a maximum in some temperature interval and a hrst order phase 
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transition, that is qnite known. Obviously, expanding the expression (^) at 
p —0 and retaining three hrst terms, one reproduces the term of Eg. (|6^ . 

This consideration convinces us that there are no discrepancies with the actual 
results of the average action method. The expression (^) gives the potential (|6^) 
predicting a hrst order phase transition. 

Our hnal remarks are on the Goldstone theorem at hnite temperature. In 
fact, at T 7^ 0 the Goldstone bosons should not inevitably be massless as it 
was argued, in particular, in Ref. [^]. Formally, the reason is that at hnite 
temperature Lorentz invariance is broken and therefore the condition = 0 
does not mean zero mass of the particle in contrast to as it should be at T = 0 
(see for details Ref.||22|). We observed in the SDA that this is the case when 
the hrst order phase transition happens. In the limit N ^ oo corresponding 
to a second order phase transition the Goldstone bosons are massless, as it is 
seen from Eg. (]T7|) . The same conclusion follows from the results of Ref. when 
a second order phase transition is assumed. Probably this problem requires a 
separate investigation by means of other methods. 

One of the authors (V.S.) thanks Institute for Theoretical Physics University 
of Leipzig for hospitality when the hnal part of this work has been done. 
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